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Bekenstein's Tensor- Vector-Scalar (TeVeS) theory has had considerable success in explaining var- 
ious phenomena without the need for dark matter. However, it is difficult to observationally discern 
the differences between TeVeS and predictions made within the A-cold dark matter concordance 
model. This implies that alternative tests are required that independently verify which theory is 
correct. For this we turn to the strong-field regime of TeVeS. fn particular, we solve the spherically 
symmetric equations of hydrostatic equilibrium for a perfect fluid with a realistic equation of state 
to build models of neutron stars in TeVeS. We show that causality within the neutron star is only 
maintained for certain cosmological values of the scalar fleld, which allows us to put constraints on 
this value independently of cosmological observations. We also discuss in detail the internal struc- 
ture of neutron stars and how each of the free parameters in the theory effects the overall size and 
mass of the neutron stars. In particular, the radii of neutron stars in TeVeS can significantly differ 
from those in General Relativity for certain values of the vector field coupling, which allows us to 
also place extra constraints on this parameter. Finally, we discuss future observations of neutron 
stars using both the electromagnetic and gravitational wave spectrums that will allow for tests of 
the appropriate theory of gravity. 
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I. INTRODUCTION 

Modified Newtonian Dynamics (MOND) [1] was de- 
veloped to describe the differences between the observed 
masses of galaxies and clusters of galaxies with the 
masses inferred from Newtonian dynamics. To this task, 
MOND has been extremely successful (for a review see [2] 
and references therein), however it is clear that MOND 
was never more than a toy model for gravitation as it is 
not a covariant theory. A relativistic version has been 
put forward by Bekenstein [3j, where the standard Ein- 
stein tensor field of General Relativity (GR) is coupled 
to a vector field as well as a scalar field, hence the theory 
is called Tensor- Vector-Scalar (TeVeS). 

In the weak acceleration limit, TeVeS reproduces 
MOND, and hence carries forth all the successes of that 
theory. In the Newtonian limit, the parametrized post- 
Newtonian (PPN) coefficients, (3 and 7, agree with the 
results of solar system tests. Moreover, on cosmologi- 
cal scales the theory has been tested against a host of 
gravitational lensing data [H , and can also be shown 
to reproduce key features of the cosmic microwave back- 
ground radiation [S] as well as galaxy distributions [7]. 
Whilst the reproduction of known observations is obvi- 
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ously critical to the success of any theory, one also desires 
predictions that differ from the standard theory, which 
here is the ACDM concordance model. To this wc turn to 
the strong-field limit in order to study both black holes 
and neutron stars. 

Black holes in TeVeS were first studied by Giannios 
[8], who solved the field equations for static, spheri- 
cally symmetric spacetimes in vacuum, hence deriving 
the Schwarzschild-TeVeS solution. He found two dis- 
tinct branches of solutions dependant on the form of 
the vector field. When the vector field is aligned with 
only the temporal direction, Giannios was able to show 
that the PPN coefficients were identical to those of a 
Schwarzschild black hole in GR. Therefore, in this regime 
the predictions of TeVeS are observationally identical to 
GR. The second class of solutions studied by Giannios 
also included a non-zero radial component to the vector 
field. Here, he was able to show that the PPN parameters 
differ from those predicted by GR. Sagi and Bekenstein 
[9J then extended this line of work and showed the so- 
lution of a static, spherically symmetric system with the 
inclusion of electromagnetic fields leads identically to the 
Reissner-Nordstrom spacetime when the vector field has 
a vanishing radial component. Moreover, Sagi & Beken- 
stein further showed that black hole thermodynamics of 
these spacetimes are identical to those in GR. To date, to 
the best of our knowledge, the equations of hydrostatic 
equilibrium in TeVeS have not been presented, and there- 
fore there is no consistent study of compact stars. 

The study of compact stars in alternative theories of 
gravity is of critical importance. Indeed in many theories 
such as f{R) gravities and scalar-tensor theories, one can 
derive observable deviations from GR that can be pre- 
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dieted with electromagnetic observations (for a review see 
[lOj). It is therefore crucial to the prosperity of Beken- 
stein's TeVeS that compact stars are studied such that 
definitive observational predictions can be made that dif- 
ferentiate the theory from its counterparts. 

In this paper we derive the TeVeS Tolman- 
Oppenheimer-Volkoff (TOV) equations of hydrostatic 
equilibrium for a spherically symmetric perfect fluid. As 
a first approximation we consider spacetimes where the 
vector field has only a non-vanishing temporal compo- 
nent. The equations are expressed in such a way that 
the deviations from both GR and scalar-tensor theories 
are obvious. Applying a realistic polytropic equation of 
state allows us to numerically solve the system of equa- 
tions, thereby deriving realistic neutron star models. In 
the strong-field limit, TeVeS has three free parameters, 
which are the strengths of the vector and scalar field cou- 
plings, K and k respectively, as well as the cosmological 
value of the scalar field tpc- We argue that the require- 
ment of causality can place restrictions on the scalar field 
such that </Jc ^ 0.001. Moreover, we use current observa- 
tions of neutron stars to constrain K <1. 

The paper is structured as follows: In section |TT] we 
provide a brief description of the main TeVeS field equa- 
tions. We do not provide a detailed derivation of these 
equations, but instead point the interested reader to ref. 
[3j . In section |III| we derive the TOV equations of hy- 
drostatic equilibrium and further discuss limits of the 
equations as the theory tends towards scalar-tensor the- 
ories and GR. In section IIVI we construct neutron star 
models, in particular discussing the equation of state in 



section IV C before going on to vary the cosmological 
value of the scalar field, vector and scalar coupling pa- 
rameters in sections |IVD[ |IV E| and |IV F| respectively. In 
section [V] we discuss the observational consequences of 
the TeVeS models as distinct from GR, and finally con- 
clude in section |VI[ Unless otherwise specified, coordi- 
nates used are x^^ = {t^r,9,(j)), where Greek indices run 
over ... 3. Antisymmetrization is denoted by square 
brackets, i.e. Aj^^j := A^^, - 
denote symmetrization, ^(^i/) 



Ai,^, and round brackets 



■A., 



II. TEVES EQUATIONS 
A. Field Equations 

TeVeS is based on three gravitational fields; an Ein- 
stein metric, (7^^, which has a well defined inverse, g^'^ , 
a timelike vector field, Ufi, and a dynamical scalar field 
ip. The physical metric, 3^,^, is defined in terms of the 
Einstein metric, the vector and scalar fields according to 



- ^-2y 



(1) 



The field equations are then derived by varying the total 
action of the theory (for details see ref. [3]) with respect 
to the Einstein tensor, the vector field, the dynamical 



scalar field and also a non-dynamical scalar field, a. By 
varying the action with respect to both scalar fields, a 
and one can show 

[^Ji{keh'^''^,^^^p) /l^V,7];5 

= kG [g'^P + (1 + U'^U^] f^p, (2) 

where a semi-colon denotes covariant differentiation with 
respect to the Einstein tensor and h^i, :— g^i, — U^Uu- 
A positive dimensionless scalar, fc, and another positive 
scalar of dimension length, £, have been introduced that 
act to couple the scalar fields to gravity. Moreover, ^{x) 
is defined with respect to a free-function J-{fJ,), which 
is not predicted by the theory. The relation between 
these two functions is given by the following differential 
equation; 

We note here that Bckenstein [3 has assumed a specific 
functional form for which led him to derive k — 0.03 
using Newtonian limits of TeVeS. As we will show, the 
results in this paper are completely independent of the 
functional form of J-'. As such, in section IV F we will 



explore the parameter space around k in an attempt to 
provide a model independent constraint on its value. 

Defining a as a "non-dynamic" field implies it has no 
kinetic term in its equation of motion. Therefore, the 
field equation on a is simply an algebraic relation; 



fcGcr^ ^ n ikfh"''^ 



V.a<^,/3) • (4) 

In this paper we only consider the strong-field limit, 
which has /i = 1 as an excellent approximation (see the 
detailed discussions in [3j|8l[9]). Therefore, we are con- 
sidering a regime where 



1 

kG' 



(5) 



Moreover, in this limit one finds that the free function 
J- does not contribute to the tensor field equations (|8|, 
and therefore it is ignored from the remainder of this 
discussion. 

In TeVeS, the timelike vector field is normalized, which 
is given by 

g'^^UMn = -1. (6) 

Note that the vector field is raised and lowered with re- 
spect to the Einstein metric. The coupling of the vector 
field to the other fields is given by the positive, dimen- 
sionless constant K. One can show that taking the limit 
as K and k go to zero implies the theory reduces to GR. 
Varying the total action with respect to the vector field 
yields 



=87rG il-e 



(g^^'-uPf^p + U^'^UPf^p) , (7) 
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where the normalization criteria ([6| has been implicitly 
used in the derivation. The use of this normalization 
implies that contraction of equation ([T]) with respect to 
the vector field gives zero identically. This implies that 
equation ([t]) is spatial, in the sense that it is orthogonal 
to the timelike vector field, U^. Moreover, it implies that 
there are just three components of equation Q, where 
the fourth component is given by the normalization con- 
dition. 

We note here that despite the extensive literature sur- 
rounding the theory of TeVeS, there is surprisingly little 
information regarding constraints on the vector coupling 
parameter, K. To the best of our knowledge, the only dis- 
cussion is from Sanders , who has phenomenologically 
argued for K < 10~^ based on solar system experiments, 
although he has not gone into detail with the calculation. 
If K is indeed this small, then the differences in neutron 
stars in TeVeS from those predicted in GR is not sig- 
nificant enough to be detectable. We further note here 



that strange things occur when K >2 (see section IV B I , 
and we therefore tend to keep the rest of our discussions 
below this value. 

Finally, variation of the total action with respect to the 
Einstein metric yields the modified Einstein equations 



Gai> ~ SttG 



T, 



©Ml'' 

(8) 



where r^,y is the physical energy-momentum tensor field 
which is defined using the physical metric and 



1 



Q/3 



AP 



(9) 



(10) 



- (1 - e-*'^) f^pU^U^ 

B. Fluid equations 



We wish to consider here a perfect fluid, which is de- 
scribed in the physical frame. Therefore, the energy mo- 
mentum tensor is 



(11) 



where u^, p and P are the fluid four- velocity, energy- 
density and isotropic pressure measured in the physical 
frame. In the next section we will choose the fluid four 
velocity and the timelike vector field to be aligned. In this 
case, Bekenstein [5] showed u 
implies 



^ — e^Uf^i, which further 



pU^U, + e-^'^Pig^,+U^U,) 



(12) 



Finally, conservation of energy-momentum is given by 

= 0, (13) 

where a vertical bar represents covariant differentiation 
with respect to the physical metric (i.e. .g^j^io- = 0) and 

^ V ■ — y av ■ 

III. HYDROSTATIC EQUILIBRIUM 

The aim of this paper is to derive the equations of 
hydrostatic equilibrium that are solutions of the TeVeS 
field equations outlined in the previous section. For this 
we shall assume a static, spherically symmetric form for 
the Einstein metric . Moreover, for this paper we shall 
choose to work in Schwarzschild coordinates, as opposed 
to isotropic coordinates which seem to be the popular 
choice in TeVeS. We do this so that the equations reduce 
to the familiar form of the TOV equations in the general 
relativistic limit. The metric is 



(14) 



where ds^ := gapdx°^dx^ is the Einstein line element and 



is the line element for the two 



3^/3 

dn"^ d6l2 + sin' 
spheres. Moreover, the spacetime being static and spher- 
ically symmetric implies all functions depend only on the 
radial coordinate. 

Analogously to choosing an ansatz for the tensor field, 
we must also select an ansatz for the vector field. In 
the cosmological limit the vector field must only have a 
timelike component such that the theory does not pick 
out a preferred direction in space. However, near strong 
gravitational fields this is not required. Instead, the only 
requirement is that the asymptotic limit of any spatial 
components must vanish. Having said that, out of math- 
ematical convenience we choose in this article to only 
consider a vector field where the spatial components van- 
ish. It is noted that this is restricting our solution space, 
and a future direction for this research is to reintroduce 
a non-zero radial component into the vector field. 

With this choice of vector field, the normalization cri- 
teria (|6| , together with the Einstein metric ( 14 1 uniquely 
determines the form of the vector field to be 



-""Z^, 0, 0, 



(15) 



In section |II A| we showed that the vector field equation 
([t]) only provides three components of the vector field, 
and the fourth is given by the normalization criteria. 
This implies that, given the form of the vector field above, 
equation Q is trivially satisfied, and indeed a simple cal- 
culation reveals this to be true. 

Utilizing the form of the vector and tensor fields along 
with the definition of the physical metric ([ij , it is trivial 
to show that the physical line element is 



ds 



^2 



-e^+^'Pdt^ + e^-^'^dr^ + e'^'^r^dn'^, 



(16) 
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where ds^ := cjapdx^dx^ . 

The conservation equations can now be expressed in 
terms of the metric coefficients. Only the radial com- 
ponent of equation (f3l is non-trivially satisfied, which 



is Euler's equation describing the pressure gradient as a 
function of the metric coefficients 



-P'=i^- + V'^(p+P), (17) 

where a prime denotes differentiation with respect to the 
radial coordinate. 

The scalar field equaiton Q can be shown to reduce 

to 

(p+iP). (18) 



-('^+0/2 



Integration of equation ( 18 ) yields 



475-7-2 



(19) 



where the "scalar mass" M^, has been defined according 
to P] 

M^{r):^4TTj (^p + 3P^ eS'^+'^^/^-^'^r^dr. (20) 

The only remaining field equations are the modified 
Einstein equations given by (|8|. For completeness, the in- 
dividual components of the Einstein equations are given 
in appendix [l] After much algebra one can show the 
modified Einstein equations reduce to 



rC -f e'^ - 1 =87rGpe^"2v^2 



Kr^ 



+ 4 V2 

1 ^SnGe'^-^'PPr^ 



k ^ 



Av' 



Kr^ 



-{2v' -2C'~v'C'r + 2v"r + v'^r) 



=8TTGPe'^-^'^r^ 



Kr^ 



, (21) 



(22) 



(23) 



We note that in GR (which is recovered by taking the 
limit K = k — 0), equation (21 1 is both independent 
of the metric coefficient, v, and is also first-order in all 
derivatives. In TeVeS however, the dependance of equa- 
tion (21 1 on v, and also the second-order terms are due 



to the presence of the vector field. These extra terms 
complicate the equations considerably, as will be seen in 
more detail in the following section. 



A. TOV equations 

We desire a form of the above equations that resemble 
the familiar TOV equations of hydrostatic equilibrium of 



GR. To this end we define a new function, m(r), accord- 
ing to 

2m 

r 



1 



(24) 



By adding equations ( 22 1 and ( 23 ) , and substituting 



the result through equation (21 ) to eliminate the second 



order derivatives in i^, one finds 



m' 1 - 



K 



2 

27rr2 



Km 
~2^ 



= AiTGe-^'^r^ (p+2KP 



Krv' 



k 



-if,- 



1 



rv 



2m\ , , 



Moreover, substituting equation (24) into equation (22 1 
yields 



,, Krv'\ SirGe-'^'^Pr'^ + 2m/r Anr ,0 
+ — ril-2m/r) +1^^ 



(26) 



Equations (17 1, (251 and (26 1, together with the equa- 



tions describing the scalar field, (19) and (20), now form 



the coupled system of equations describing hydrostatic 
equilibrium in TeVeS. We will show in the next section 
that these reduce to the standard TOV equations of both 
scalar-tensor theories and GR in the appropriate limits. 
Equation ( 26 1 is expressed in an illustrative form such 



that the differences between GR and TeVeS are obvious 
(see the discussion in section IIIBl. However, we note 



that equation ( 26 1 is not expressed in a convenient form 
due to the quadratic term in v'. We can therefore al- 
gebraically solve this equation for v' , such that it is ex- 
pressed explicitly; 



Kr 



±<l + K 



AnGe-'^'fi Pr^ + m/r 2nr'^ 



r (1 — 2m/r) 



k 



1/2 



(27) 



One can see that there are two solutions to the above 
equation, however we note that the negative sign is not 
physical. In this case it is trivial to see that the limit 
K ^ does not exist. We therefore exclude this solution 
from the remainder of the discussion. 

Before moving on, there exists another convenient 
equation that will subsequently be used that is derived 



from substituting equations ( 26 1 through ( 25 ) to yield 



m 
r 



V r 



2m 
r 



l6nGe-^'^Pr^ 



AnGe-^'^r^ 



(p + 3P) 



K 
~2 

(28) 



From the grouping of the terms on the left hand side of 
this equation, it can be used to learn much about the 
range of the vector field coupling constant, K. We anal- 



yse this in more detail in section IV B 
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B. Limits of equations 



A. Exterior region 



1. Scalar-Tensor theory 

TeVeS couples the vector field to the other fields 
through the positive constant K. Setting this to zero 
is equivalent to ignoring the eff'ects of the vector field in 
the equations. This is best seen through the definition of 
the vector field action, which is directly proportional to 
K (see eqn. (26) of Bekenstein [3]). 

With K = Q, the quadratic term in v' in equation (26) 
vanishes. Moreover 



equation ( |25[ ) loses the explicit cou- 
pling between the function m and the metric coefficient 



m! = AirGpe 



1 - 



2m 



STrGe-^'fpr^ + 2m/r Anr 



7' (1 — 2m/r) 



k 



-if- 



(29) 
(30) 



These equations are now the same as the well known 
TOV equations for scalar-tensor gravity [T^ I13j . 



2. General Relativity 

The GR limit is reached by also letting k = 0. We 
note from equation (19) that (p' is proportional to k^, 



implying ip'/k oc k. In this limit, equation (19) can be 
integrated to give ip — ipc which is constant. Then 



m =47rGe '^'^'^ / /5r^dr, 
Jo 

, ^TiGe-'^f-Pr'^ + 2ml7 
V =- 



r (1 — 2m/r) 



-2P' 
'p+P 



(31) 
(32) 

(33) 



which are the standard TOV equations of GR with a con- 
stant, background scalar field, (pc- One can always choose 
an appropriate coordinate transformation such that this 
scalar field vanishes, implying the equations reduce to 
the standard TOV equations in GR. Moreover, equation 



(31 1 is a usual definition for mass, and we therefore refer 
to the function m as a mass function, even in the TeVeS 
regime where it is not clear whether this truly represents 
any physical mass of the system. 



IV. NEUTRON STAR MODEL 

In this section we concentrate on realistic neutron star 
models based on the full set of TOV equations in the 
TeVeS theory derived in the previous section. We begin 
by briefly looking at the spacetime external to the neu- 
tron star in order to derive observable quantities, such as 
the ADM mass of the neutron star. 



The boundary of a neutron star, denoted r^, is defined 
as the point at which the pressure vanishes, P — 0. For 
any r > rs (for any reasonable equation of state; see 
below) , this implies that the density is also zero, and the 
spacetime is therefore given by the Schwarzschild-TeVeS 
solution. We note that an analytic solution with P — 
p = in Schwarzschild coordinates is possible, however 
its implicit nature provides no new information to the 
present discussion, and hence we omit its inclusion. An 
analytic vacuum solution in a nicer form is known in 
isotropic coordinates, and was first derived by Giannios 
[8J. For this discussion we simply require the asymptotic 
behaviour of the vacuum equations far from the source, 
such that we can gain an understanding of the function 
m and the scalar field in the asymptotic limit, and also 
derive the ADM mass of the system. 

For deriving observable quantities in TeVeS we wish to 
use the physical metric. However, the physical metric. 



as expressed in equation (16), is not in an asymptoti 



cally fiat form. Therefore, to derive the ADM mass for 
the system being considered herein, we require a coordi- 
nate transformation that brings equation (16) into such 



a form. Bearing in mind that the cosmological value of 
the scalar field is defined as the limit of the scalar field 
as r — > oo, we perform a coordinate transformation given 
hy t = te'^" and f = which yields 



ds 



^2 



(34) 



Setting p = P = in the TOV equations gives the 
Schwarzschild- TeVeS solution in asymptotically flat co- 
ordinates. By expandin g m and v in infinite series', one 



can integrate equation ( 19 
plies is constant) whic 



(noting that p — P 
ti gives 



kGM, 



O 



1 



im- 



(35) 



which is in agreement with [S]. Further substituting the 
series expansions through the remaining equations im- 
plies the asymptotic form of the Schwarzshild-TeVeS so- 
lution can be shown to be 



9 it = -1 
5f f = 1 + 



2e~f- 

1 — 

r 

2e-'f- 



kGM^ 

kGM^ 
An 



O 



1 



O 



Therefore, the ADM mass of the system is [14] 



Mabm = rri, 



kGM^ 
An 



(36) 
(37) 

(38) 



where moo is the mass function evaluated at radial infin- 
ity. As expected, this implies that the scalar field mass 
contributes to the entire mass of the spacetime. In GR, 
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one finds that the mass function outside the star is a 
constant, and therefore the ADM mass equals the mass 
function for all r > r^. However, in the TeVeS case, the 
non-zero contribution of the scalar field implies that the 
mass function still varies outside the star (see figure |4]) , 
and that this is not equivalent to the ADM mass of the 
spacetime. 



B. The coupling parameters 

We can utilize the equations derived herein to put sig- 
nificant constraints on the parameter K. Firstly, by look- 



ing at equation ( 28 1 , one can see that the structure of the 
equations significantly changes either side oi K = 2. In- 
deed, for K = 2\t can be shown that an equation of state 
is forced upon the system, namely p — —3P. Moreover, 
analysing the remaining equations shows that the system 
has become underdetermined due to a redundancy in the 
equations. 

We can analytically determine that values oi K > 2 
are also unphysical. This can be seen by performing a 
series expansion of the hydrostatic equilibrium equations 
near r — 0. After much work, one can show 



P^Po 

z • ■ 

m ^msr^ + O (r^) , 
where the expansion coefficients are given by 



^2 

P2 
TO3 



kG 

"3" 

167rG 



{2ifi2 + 1^2), 



3(2-K) 

Po + pQ 
2 

AttG 

''3{2-Kf 



— 2<po 



2po + 3KPv. 



(39) 
(40) 

(41) 
(42) 

(43) 
(44) 

(45) 
(46) 



It is trivial to see that K > 2, and any reasonable equa- 
tion of state (i.e. such that p+3P > 0), implies < 0, 
and hence m < for small radii. It is not entirely clear 
whether a negative mass function itself is unphysical. 
However, it can be shown that for models with K > 2, the 
pressure diverges from the stellar centre outward. This 
implies a stellar model is not possible as the pressure 
never vanishes, and hence the surface of the star is not 
defined. Therefore, we can restrict the values of the vec- 
tor coupling parameter to < < 2. This is consistent 
with Sagi and Bekenstein [S], who showed that physical 
black hole solutions are only valid for X < 2. 



C. Equation of State 

As mentioned, the TOV system of equations are closed 
with an equation of state relating the pressure and the 
density. There has been considerable interest over recent 
years on the neutron star equation of state (see (TS] and 
references therein). If GR is the correct theory of grav- 
ity, then information can be gained on the appropriate 
equation of state from the masses and radii of neutron 
stars. However, when one considers alternative theories 
of gravity there is a degeneracy between the equation of 
state and the free parameters of the theory. Indeed this 
is true in TeVeS as we will show below. 

In order to concentrate our attention on the structural 
dependence of neutron stars on the free-parameters in 
TeVeS, we choose to only work with a single equation of 
state in this paper. It is noted however that this implies 
we can not stringently constrain these parameters due to 
the intrinsic uncertainties associated with our choice of 
equation. We use a polytropic equation of state given by 



p 


=JCnQmb ( — ) , 


(47) 


p 


P 


(48) 


rrib 


=1.66 X lO"^'* g, 


(49) 


no 


=0.1 fm"^ 


(50) 



where T = 2.46 and /C = 0.00936. Here, the parameters 
r and K, have been fitted to the tabulated data for the 
realistic equation of state, known as equation of state A 
given in Arnett and Bowers |16j . This equation of state 
has been extensively used in previous studies of neutron 
stars in GR and scalar-tensor theories (see for example 
Damour and Esposito-Farese [12^, Sotani and Kokkotas 

The interior conditions for the neutron stars are found 
using equations ( 39p6 1. Using these conditions we can 
determine stellar models for various parameters associ- 
ated with the theory by numerically solving the TeVeS- 
TOV equations derived in the previous sections with the 
aforementioned equation of state. Numerically, we con- 
sider radial infinity to be r = 100 km, and have verified 
that the results do not change when we expand our nu- 
merical box to r = 500 km. The structure of neutron 
stars depends on the various parameters in the TeVeS 
theory, for example the vector and scalar field coupling 
parameters, K and k respectively, and also the cosmo- 
logical value of the scalar field, ipc- We therefore spend 
a section on the variation of each of these parameters. 



D. Cosmological Value of the Scalar Field 

Cosmological considerations imply the asymptotic 
value of the scalar field is restricted to < ipc 1 
[3]. However, using neutron stars we can place a further 
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restriction on the lower bound for this value, as shown 
presently. 

Bekenstein showed that TeVeS allows for superluminal 
propagation of tensor, vector and scalar perturbations 
when (f < (see section VIII of [3^). This was shown 
using perturbations of the various fields in the physical 
frame, and was independent of the matter content of the 
model. Therefore, for neutron star models we demand 
that the scalar field be everywhere greater than or equal 
to zero. This can be used to constrain the cosmologi- 
cal value of the scalar field. Figure [T] shows the scalar 
field as a function of the radius for various values of the 
vector coupling parameter. Here, each line corresponds 
to the stellar model with the maximum mass (all figures 
henceforth that are functions of the radius also corre- 
spond to the models with the maximum mass). In fig. 
[l] fc = 0.03, the top series of lines have ipc = 0.003 and 
the bottom are for ipc = 0.001. One can see that for 
(fie = 0.001, models where K < 1 are unphysical as the 
interior of the star has (p < 0. However for larger values 
of ifc, we find that all models are physical, independent 
of the vector coupling parameter, K. In section [IV E| we 
show that K <1 based on considerations of the mass of 
neutron stars. Therefore, from the demand for sublumi- 
nal propagation of perturbations, we can conclude that 
ipc ^ 0.001. We note that this is a conservative estimate, 
and a true minimum value may be higher than this due to 
variations in the equation of state. Moreover, if one can 
further constrain the value of K to be smaller than unity, 
then the minimum cosmological value of the scalar field 
will also rise accordingly. We also mention the caveat 
that we have assumed only a temporal component of the 
vector field which may alter these results. 
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FIG. 1: Scalar field as a function of radius with fixed k = 
0.03. Here, the top four lines have = 0.003 and the bottom 
four lines have = 0.001. We exclude all models where (^3 < 
for any value of r as this implies violation of causality. 



Whilst in the above we assumed that fc = 0.03 accord- 
ing to Bekenstein [3], it is worth exploring the variation 
of these values as a function of the scalar field coupling. 
Figure [2] again shows the scalar field as a function of the 
radius, however this time we hold K fixed at 0.1 and vary 



the value of fc. One can see that when ip^ = 0.001, only 
very small values of the scalar coupling permit models 
where > for all r. However, when p^. — 0.003, mod- 
els for fc < 0.04 are accepted. 




FIG. 2: Scalar field as a function of radius with fixed K — 0.1. 
The top three lines have ipc ~ 0.003 and the bottom three 
lines are for ipc = 0.001. 



An important quantity that will be considered 
throughout this article is the mass-radius relation, where 
we are using the observable ADM mass, Madm, as a 
function of the physical radius, R := e^'^'^''^^^. Figure [s] 
shows the dependence of ipc on the structure of neutron 
stars for two different values of K. One can see here that 
the dependence on the cosmological scalar field on the 
size of neutron stars is completely negligble. This implies 
that, whilst ipc may evolve with cosmological time, the 
overall structure of neutron stars throughout this evolu- 
tion will not significantly change. 
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ADM mass as a function of the physical radius, 



Here, k 



0.03 and we show two different 
values of the vector coupling parameter, K = 0.5 and 1.0 for 
different cosmological values of the scalar field. 
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E. Vector Field Coupling 

The vector field coupling is possibly the least restricted 
parameter in the TeVeS theory. To the best of our 
knowledge, the only discussion of this parameter is given 
by Sanders [TT], who postulated on phenomenological 
grounds that K < 10"''. In this article we will pro- 
vide constraints on the vector field of K < 1 based solely 
on observations of neutron stars. It is difficult to fur- 
ther constrain the parameter using this method due to 
uncertainties in the equation of state. 

Figures |4] and |5] show the structure of the mass func- 
tion, m, and the metric coefficient, v, in the GR case 
as well as for different values of K, with k = 0.03 and 
(fie — 0.003. The stronger the vector coupling (i.e. higher 
K) implies the mass of the neutron star is decreasing. We 
note that the mass function has a local maximum which 
is different from the asymptotic value as one would ex- 
pect. This is because the scalar field acts to decrease 
the overall mass of the star exterior to the stellar surface 
(this was also discussed briefly in Giannios f^). 
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FIG. 4: Mass function, m, as a function of radius with k 
0.03 and ipc = 0.003. 



In figure [6] we show the dependence of the measurable 
ADM mass on the central density of the star, po := p{r = 
0), for varying values of the vector coupling parameter, 
K. Figure [7] shows the ADM mass-radius relation of the 
star for the same models as in figure [6j In both these 
figures we use k — 0.03 and fc = 0.003. One can see 
from these two figures that for small values of K, i.e. 
K < 0.05, the deviations from GR are small, and indeed 
these effects are difficult to distinguish between varia- 
tions around the equation of state and different values of 
K. However, for larger values oi K > 0.05, the devia- 
tions from GR begin to become considerable. Indeed, in 
these models one can not create neutron stars with mass 
greater than approximately 1.3 Mq. When the value of 
K is greater than unity, neutron stars begin to look ex- 
tremely different than those predicted in GR. We find 
that, independent of the central density, one can not cre- 
ate a neutron star with mass greater than approximately 




FIG. 5: Metric coefficient, v, as a function of radius with 
k = 0.03 and (fie = 0.003. 



one solar mass. Moreover, the radius in these models is 
still less than approximately seven kilometres. 



= 0.03 0.1, 
ip^ = 0.003 




FIG. 6: Central density, po :— p{r = 0) as a function of the 
ADM mass, Madm for varying values of the vector coupling 
constant K. Here, the scalar coupling constant is held con- 
stant at = 0.03 and the cosmological value of the scalar 
field is ifc = 0.003. Variations for the ADM mass are small 
for small values of the vector coupling constant K, however 
they grow as K tends towards 2. 



F. Scalar Field Coupling 

The scalar field coupling has been constrained by 
Bekenstein [3] using planetary motions in the outer solar 
system to be fc w 0.03. However, we note that this used 
a specific functional form of the free function J- , whereas 
in this paper our results have been shown to be inde- 
pendent of this function. Therefore, in this section we 
explore the parameter space of k in an attempt to place 
model independent constraints on its value. 

Figure [8] shows the mass-radius relation for different 
values of the scalar field coupling parameter. Here, (pc = 
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FIG. 7: Mass-radius relation for various values of the vector 
coupling parameter K. Again, k — 0.03 and tpc = 0.003. For 
small values of K there is little difference in the mass-radius 
relation from the predictions from GR. However there are 
significant differences for larger values of K, where the mass 
of a neutron star with the equation of state being considered 



here must be smaller than Madm ~ 1 M( 



0- 



0.003 and K = 0.5 and 1.0. One can see that that value 
of the scalar coupling parameter needs to be significantly 
altered to change the size of the neutron stars. It is 
interesting to note that increasing the value of k creates 
smaller neutron stars. As mentioned, some models with 
larger values of k may be unphysical due to the causality 
problem discussed in section [IV D[ However, considering 
fig. [3] we can still conclude that the dependence of the 
scalar field coupling, k, on the structure of neutron stars 
is not so strong. 




6 7 8 

R [km] 



9 10 



FIG. 8: Mass-radius relation for variations of the scalar field 
coupling parameter. Here, ipc ~ 0.003 and the vector coupling 
parameter is K = 0.5 and K = 1.0. 



V. OBSERVATIONAL CONSEQUENCES 

Thus far, we have investigated the structure of neu- 
tron stars in TeVeS as distinct from GR. It is worthwhile 
spending some time establishing how this can be seen 
observationally, and therefore provide independent veri- 
fications or falsifications of TeVeS. We have seen in the 
above sections that the biggest variations in neutron star 
structure are due to the vector field coupling parameter, 
K. Therefore, throughout the remainder of this section 
we choose to leave k = 0.03 and (pc — 0.003, and only 
vary the values of K. It is worth noting at this point 
that if the discussion of Sanders [TTj is correct (i.e. that 
K < 10^^), then the differences in neutron star predic- 
tions by GR and by TeVeS are observationally indistin- 
guishable. 

Figures [6] and [7] suggest that the mass and radius of 
neutron stars in TeVeS are smaller than those in GR. 
That is, a higher value of the vector coupling constant, 
K, implies smaller masses and radii for neutron stars. 
The most accurate measurements of neutron star masses 
are from timing observations of binary pulsars, although 
significant progress is also being made using binaries with 
white dwarf companions (see Lattimer and Prakash |15j 
and references therein) . Conservatively, one can say that 
neutron stars exist with masses greater than ~ 1.5Mq, 
although the true value may be significantly greater. 
Taking this into account, together with the uncertainties 
in the equation of state that we have been using, suggests 
we can rule out the parameter space K > 1, although we 
note that this too is probably a conservative estimate. 
Moreover, as TeVeS systematically predicts lower masses 
for a fixed equation of state, if a set of equations of state 
are excluded by the measurement of a high mass neu- 
tron star in GR, then these equations of state are also 
excluded in TeVeS. For example, if the mass of a neu- 
tron star is measured at M > 1.7Mq, then our adopted 
equation of state is ruled out in both GR and TeVeS. 

Another observable quantity is the redshift of atomic 
lines emanating from the surface of the neutron stars. 
The redshift is given by 



(51) 



where gtt,oo and gtt,s are the values of the temporal com- 
ponent of the physical metric at radial infinity and the 
stellar surface respectively. With this redshift factor z, 
the photon energy emitted from the stellar surface. Eg, 
can be connected to the energy measured at infinity, Eao , 
as = [1/(1 + z)]E, [ig. 

Figure [9] shows the redshift at radial infinity of a pho- 
ton emitted from the surface of the neutron star as a func- 
tion of the mass for various values of K with fc = 0.03 and 
(fie = 0.003. Again, for small values of the vector coupling 
parameter the deviations from GR are small. For larger 
values however, one finds a large redshift from extremely 
small values of the mass. That is, for neutron stars of 
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Madm ^ 1 Mq, one can still have redshift z 0.4 
1 
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FIG. 9: Redshift as a function of the ADM mass for various 
values of K, with k — 0.03 and ipc — 0.003. 

An exciting future probe of neutron stars and black 
holes is through gravitational wave astronomy. Neu- 
tron stars in metric theories of gravity are known to 
exhibit pulsations of the spacetime, known as w-modes 
[19j . These w-modes have frequencies, where /^i? cx 
Madm/R (see Kokkotas and Schmidt j20j for a discus- 
sion in GR and Sotani and Kokkotas [H] for scalar-tensor 
theory) . Madm / R is known as the compactness, and we 
plot this against the ADM mass for different values of 
the vector coupling, K, with k = 0.03 and (pc — 0.003 in 



figure 10 If future gravitational wave detectors are able 
to reach the required sensitivity to accurately determine 
the compactness of neutron stars through w-mode oscil- 
lations, then these results coupled to the mass estimates 
from pulsar timing observations will be able to provide 
an independent constraint on the vector coupling con- 
stant, K. Moreover, it has been shown that the relation- 
ship between the frequency of the w-mode and the com- 
pactness is almost independent of the equation of state 
[20|, I21j . This implies that gravitational waves will pro- 
vide a unique probe of the appropriate theory of gravity 
which, unlike most other observations of neutron stars, 
does not depend on the equation of state. 



VI. CONCLUSION 

We have solved here the equations of hydrostatic equi- 
librium for a spherically symmetric, static distribution 
of a perfect fluid in Bekenstein's Tensor- Vector-Scalar 
(TeVeS) theory [3 . Further imposing a specific poly- 
tropic equation of state has allowed us to analyse the 
differences in neutron star structure between those ap- 
pearing in General Relativity and those in TeVeS. We 
have shown that the form of neutron stars in TeVeS de- 
pends most strongly on the coupling parameter associ- 
ated with the vector field, and used this to show that 
K < 1. Constraining this parameter any further is diffi- 
cult with current observations given the uncertainties in 
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FIG. 10: Compactness, Madm/R, as a function of mass for 
various values of K, with k = 0.03 and fc = 0.003. 



the equation of state. We have further shown that the 
cosmological value of the scalar field, (pc, is constrained 
in this theory by demanding that causality not be vio- 
lated in the interior of neutron stars. In some ways, this 
is a counter-intuitive result, as the global, cosmological 
value is determined by local measurements. This led to 
the conservative constraint that (pc > 0.003, where the 
uncertainty is again due to the equation of state as well 
as the choice of vector field only having a temporal com- 
ponent. 

Using the results derived herein we have looked at ob- 
servable quantities such as the mass-radius relation, sur- 
face redshift and compactness which may allow future, 
independent verifications or falsifications of the theory. 
If General Relativity is the correct theory of gravity, then 
these observations are sufficient to determine the equa- 
tion of state of neutron stars. However, in alternative 
models of gravity there exists a degeneracy between the 
equation of state and the specific free parameters in the 
theory. This implies that independent establishments of 
the neutron star equation of state are required in order to 
stringently constrain the free parameters of any theory. 

This work has provided a first step in the study of 
neutron stars in TeVeS. The first and most obvious ex- 
tension is to allow for a non-zero radial component of the 
vector field. To do this, one simply requires that the ra- 
dial component vanish at spacelike infinity such that the 
theory does not pick out a preferred direction in space. 
In this way, the vector field equation (|7| will no longer 
be trivially satisfied. Instead, this will provide one equa- 
tion, together with the normalization criteria which will 
allow for the complete determination of the vector field in 
terms of the metric coefficients. As was shown by Gian- 
nios [8] , the introduction of a non-zero radial component 
can considerably alter the structure and geometry of the 
black hole solutions, and we anticipate this to also be the 
case for neutron star solutions. 

With basic neutron star models in place, one can begin 
to study more complicated aspects of the observations in 
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order to further constrain TeVeS. For example, recent where all other components vanish. The tensor r^^, as 

observations of quasi-periodic oscillations in the giant defined in ([9| is 

flares associated with soft gamma repeaters are believed 

to be associated with the global oscillations of strongly 

magnetized neutron stars (for example see [21] ) • 

Perturbations of neutron stars in General Relativity can 

successfully account for all of the frequencies associated 

with the quasi-periodic oscillations [25 . Therefore, by e"^'~- 

performing a similar analysis on neutron stars in TeVeS tu ~~2j7q^ ' (^'^) 

one should theoretically be able to further constrain the 

various parameters associated with the theory. This work ^2kG'^''^ ' (^^) 

is obviously non-trivial as one first needs to derive the 2 

coupled system of equations governing magnetic fields in -j-^g ~ !Lg-C(p'2 _ ^^^^ (A9) 

TeVeS, at least to a first-order approximation. 2fcG sin 9 



APPENDIX 1: COMPONENTS OF EINSTEIN 
EQUATIONS 



Here we show the various components of the modified 
Einstein field equations Q which lead to the derivations 
of equations ( 2Tp3 1. Firstly, the Einstein tensor, G^, 
which is derived with respect to the Einstein metric, ( 14 1, 
has four non-zero components 



Finally, defining t^, f^, + (l - e-^^)U''f^i^^U,), it 



is found 



G,, ^^(yV-e'^ + l), 



(Al) 
(A2) 



Gee ^-^E-^ - 2C' - i^'C'r + 2u" + u'^r) . (A3) 



The tensor S^ji/, as defined in equation (10) is given in 
component form as 



K 



/2 



inG (1 - e-**^) 



K 



Bee =^ie-W = 



(A4) 
(A5) 

(A6) 



Tit = (2e-2'^ - e^'^) pe^ 
fee =e-2^r2p : 



1 



sin^ei 
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